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Abstract 

We use some properties of orthogonal polynomials to provide a class of upper/lower vari- 
ance bounds for a function g{X) of an absolutely continuous random variable X, in terms 
of the derivatives of g up to some order The new bounds are better than the existing ones. 
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1 Introduction 

Let Z be a standard normal random variable and ^ : II — )• R any absolutely continuous function 
with derivative g'. Chernoff, [8], using Hermite polynomials, proved that (see also the previous 
papers by Nash, [13] and Brascamp and Lieb [3]) 

Ym-g{Z)<E{g'{Z))\ 

provided that E(^'(Z)) < oo, where the equality holds if and only if g is a polynomial of 
degree at most one - a linear function. This inequality plays an important role in the isoperi- 
metric problem and has been extended and generalized by several authors; see, e.g., [7, 17, 5, 
15, 10, 14] and references therein. 

Definition 1.1 (Integrated Pearson Family) Let X be a random variable with density f and 
finite mean /i = EX. We say thatX (or its density) belongs to the integrated Pearson family (or 
integrated Pearson system) if there exists a quadratic polynomial q{x) = 5x^ + ^x+y (with 
5,j3,7G R, |5| + |j8| + I/I > Oj such that 

r {pL-t)f{t)At=q{x)f{x) forall ;c G R. (1.1) 

J — oo 

This fact will be denoted by X or f ^lP{iJ.; q) or, more explicitly, X or f ^W{lJL;8 ,^ ,'Y). 

Let a random variable X ^ W{pL;q) and consider a suitable function g. Johnson, [11], 
established Poincare-type upper/lower bounds for the variance of ^(X) of the form 

(-l)"(Varg(X)-S„)>0, where S,, = Yi-lf-'^^^^^^^-^^^ . (L2) 
In particular, for the normal see [10, 15] and for the gamma see [15]. 
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Afendras et al., [1], using Bessel's inequality, showed that 

Varg X > \ , ^ ^ ,f , ^ ^ ; (1.3) 

for the case n = \ see [4]. 

Afendras and Papadatos, [2], showed that, under appropriate conditions, the following two 
forms of Chernoff-type upper bounds of the variance of are valid: 

ti/!E^'(^)n;L^-'i(i-75) ("+i)!n5i;Hi-75) 

EV(X)g(')(X) Eg"(X)(g(")(X))' 

The equality in (1.4) holds when g is a polynomial of degree at most « + 1 and in (1.5) when 
g is a polynomial of degree at most n. The bound (1.5) for beta distributions is also given by 
[18]. 

In Section 2 we provide a new class of upper/lower bounds for the variance of g(X). They 
can be called as "Poincare-type" of order n and with point balance m. They hold for a subfamily 
of Pearson distributions. In particular, the bound for N{p.,G^) distribution, namely 

^.-fe) = E^EV(X) + E(-l)-^Efe<..(X))^ 

(for {x)k see Definition 2.1), satisfies the inequality 

(-l)"(Varg(X)-5„,,„(g)) >0, 

where the equality holds if and only if g is a polynomial of degree at most m + n. 

For n fixed. Section 3 investigates the bounds S,„^n{g) as m increases. It is shown that the 
bound S,„+i,,j(g) is better than S,n,,i{g), i-c, 

I Var g(X ) - 1 {g)\<\ Yai-g{X ) - (g) \ . 

Also, for any suitable function g, S„,.„(g) — Varg(X) as m — oo. 



2 Universal variance bounds 

This Section presents a wide class of variance bounds. First we give the following definitions 
that will be used in the sequel. 

Definition 2.1 For G K, and S IN define: 

(a) {x)k = x{x — 1 ) • • • (x — ^ + 1 ), with {x)q = 1. 

(b) [x]k=x{x+l)---{x + k-l), with [;c]o = 1. 
Note that [x]^ = {-l)''{-x)k = {x + k- l)^. 
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Definition 2.2 (cf. [2]) Assume that X ~ IP(jU;^) and denote by q{x) = 5x^ + px + y its 
quadratic polynomial. Let J{X) = (a; ft)) be the support ofX and fix the non-negative integers 
m,nsuch that E|Xp^ is finite, where £ = max{m,n}. We shall denote by J^'"'"{X) the class of 
Borel functions g : (a, ftj) — t- B, satisfying the following properties: 

Hi : g G C^^' (a, (o) and the function g^''^^^\x) := '^^ffl is absolutely continuous in [a, (o) 
with a.s. derivative g^^\ 

H2 : ^q"{X){g^"\X)f < 00 and ^q'{X)\g^'\X)\ < 00 for i = n+ \, ... ,m, if any. 
Note that for m = n = 0, the property Hi does not impose any restriction on g, and 

jf'O.O(X) =l2(E,X) := {g: (a,ftj) ^Rsuch thatVarg(X) < 00}. 
Also, it is obvious that J^^'" = J^^'" = ■■■ = «'". 

Furthermore, we shall denote by " "(X) and J^°°{X) = J^'"'°°{X) [m is arbitrary because 
in this case this index is insignificant] the classes of functions Jif°°'"{X) := nJ^^Q^"''"(X) and 
:= n~^oJr~'"(X); i.e, 

^~'«(X) = {ge C"(/) : ^q"{X){g'-"\x)y- < 00 and ^q'{X)\g^'\x)\ < 00 V/ > n}, 
^"(X) = {ge C~(/) : Eq"{X){g'^"\x)f < 00 V« G M}. 

From (A.S) and from EV(^)|gW(X)| < Eg'(X) • Eg'(X)(gW(X))^ < 00, / = 0, 1, . . . 
we conclude that the (finite or infinite) sequence J^"''"(X) is decreasing in m and in n. In 
particular, if all moments exist then 

l2(r,x) = je^^^ix) 

lU 

Jfi'°(X) D jri'i(x) 

lU lU 

jf2'0(x) D jr2'i(x) D jr2'2(x) 

lU lU lU 

lU lU lU 

jf°°'°(x) D jr~i(x) D jr~2(x) D ••• D ^"(x) 

LetX ~ IP(jU;5,j3,7) with 5 < 0. Also, consider two (fixed) non-negative integers m,n, 
with « 7^ 0, and a function g G J^'"'"(X). According to Parseval's identity we have that 

00 

Varg(X) = £cf, (2.1) 

k=i 

where Ck = Eg(X)<p<:(X) are the Fourier coefficients of g with respect to the corresponding (to 
X) orthonormal polynomial system {(Pk}'k=o- 

For each / = 1,2, ... ,« the function g^'^ G J^'"^''"^'(X,) and, from Parceval's identity again, 
E(g('^(X,-))^ = ZT=Q where = Eg(')(X,-)<pi,,-(X,-) ai-e the Fourier coefficients of g(') 
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with respect to the orthonormal polynomial system {(Pk.i}T=o corresponding to X,- ~ /,■ oc g'/; 
see Appendix A. Using (A.4) we have that 

=» k+i-2 . 

WiX){g('\x)f = '£Uk)i n i^-jS))cl i = l,2,...,n, (2.2) 

k=i j=k-\ 

where each coefficient of is positive. Also, from (A.2), 

E^'(Z)g«(X) = (i\Eq'{X) if (1-75))'^%, /= l,2,...,m. (2.3) 

Let A„ = ^271, . . . , A„7,)' € m". According to Tonelli's Theorem we have that 

v=imi\hn{k)iU%z\{i -j8)\ci=ru \hn\i:T=i miU)ti-Mi-jS)]ci=ru \kn\ x 

using Fubini's Theorem, 

n CO mm{k,n} k+i-2 

£A,;„E^'(X)(g«(X))'=£p,.,4, where p,.„ = £ n(l-^'^)- (2.4) 

;=1 /t=l 1=1 y=<:-l 

We seek for a vector Am,„ such that Pm+i.n = Pm+2,n = ■■■ = Pm+n,n = 1- From (2.2) we obtain 
the system of equations 

^m,n ' ^m,n — In; (2.5) 

where the matrix A„,,„ G M"^" has (r,c)-element which is given by 

and the vector 1„ = (1, 1, . . . , 1)' G E". 

The above system has the unique solution, see Appendix B, 

= (-i)'-i c;)/ [(m + «),-n;L+;r ^1 - js)] , / = 1,2, . . . (2.6) 

From (2.4) and (2.6), using the hypergeometric series (C.l), we have that Pk-m,n = 1 — + 
« - ^).n;t+:::rHl - 7-5)/ [(m + «)„n7ir Hi - 7-5)] . Equivalently, 



Pk;m,n — * 



"";!^;rf,/..:" , i<k<m 



Thus, 



(m+n-^)„n;L+';:r'(i-;g) 
(M+«)„n"Sr' (1-7-5) 

1 , m<k<m + n, 

i+(-i)- ":":';-gsC".f' - ">'"+"■ 

(")+«)« n"=„" (i-jo) 



= Varg(X) - y ^^^C + —4 (2-7) 

y (-1)"-!^ ^cf- 
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The main result of this paper is contained in the following theorem. 

Theorem 2.1 Let X ~ IP(/i; 5, j8, 7) with 5 < 0. Fixed two non-negative integers m,n [with 
n 7^ 07 and a function g € Jif"''"{X). Consider the quantity 



SmAg) = Ifl,-EV(^k^'H^) + I(-i)'"'W(^)(5*'^(^))', (2.8) 

where 



(T)n:;s:;TrHi-75) 

a,- — 



(m + n),E9'(^) (njL-' 1 (1 - j5)) njir ' (1 - js) ' 



{'n+n)tnj:,r{^-jS) 

are strictly positive constants ( depending only onm,n and X ) and the empty sums ( when m = 0) 
are treated as zero. Then the following inequality holds: 

(-l)"(Varg(X)-5,„,„(g)) >0, 

and where S,n^n{g) becomes equal to Varg^X) if and only ifg is a polynomial of degree at most 
m + n. 

Proof From (2.8), via (2.3) and (2.7), we obtain that (-1)" (Varg(X) - 5,„,„(g)) = 
Rm.n{g)^ where the residual 

R,n,n{8)= L ^tMx)ct:= L fa+») n"'+"-ifi-/g) ^ ^ ^ 

k>m+n k>m+n + ^)n\.ij=,n U ](^) 

is non-negative and equals to zero if and only if Cyi = for all k > m + n, i.e., the function 
g : J{X) — > R is a polynomial of degree at most m + n. □ 

Remark 2.1 (a) For fixed n and for any function g G J^'^-"{X), where M can be finite or 
infinite, the variance bounds {Sm.n{s)}m=o °f same type, i.e. upper bound when n is 
odd and lower bound when n is even. 

(b) The bounds {S,n,n{g)}"„=o require the same conditions on the function g, i.e., g E J^"'"{X). 
Remark 2.2 (a) The bound Si,i (g) is the bound (str) of (1-4). 

(b) The bounds So.„(g) are the bounds S„ which are given by (1.2). Also, for the special case 
m = 0,n = \ observe that 5o.i (g) = 5i = 5i (^gak) ' see (1.5). 

(c) The results shown in Theorem 2. 1 apply to the special case where n = (note that the 
second sum is empty and is treated as zero). In this case the lower bound Sm.aig) reduced to 
the one given by (1.3). 

Remark 2.3 Regarding the conditions on the function g of Theorem 2.1 we note that g G 

j^max{m,n},n-i(^X) \ niax{m,n},« (j,^ ) implies that the bouud 5„,,„(g) is trivial, i.e., +00 when n 
is odd and —00 when n is even. 
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Table 2.1: Specific form of the bounds 5',„,„(g) for normal, gamma and beta distributions. 

distribution parameters / J(X) q(x) 
bounds S,„^(g) 

normalfix.a^) ix e R, > i_e-(f-M)-/2(T- jr, (j^ 



gamma(a, 0) a,0 >0 rm^^^^'' (0,°=) 



gli^iE¥'-)(^)+LLi(-i)''-'(SJ^(g"'H^)) 


r(a) 

g'li(^;i^'EW'H^)+L;Li(-i)'-'(JjprEx^^ 

beta(a,iS) «'/^ > (1 -^)''-' (0,1) x(l -x)/(a + j8) 

ym (7)|«+/i+».+.1„,[«+/i|,,- -nr2y,7i yy-.f-Vy'i 

■^1=1 (m+«),|a|i[/5|,|a+/S+i-l];[a+j3+m]„ '-^ ^'^ 

+ILi(-i)^-' (,„+„il,+„,, iEr(i-^)'(g»W)' 



Now, we seek for upper bounds of the non-negative residual /?,„,„(§). 

Proposition 2.1 Assume the conditions of Theorem 2.1 and, further, suppose thatgG J^^'^{X) 
for some T G {«,..., wi + ?i + 1 }. Then the residual R,„,n (g), given by (2.9), is bounded above 
by 

UrEq'{X){g'^'\X)f, T = n,n+l,...,T, (2.10) 
wlwre Ur=u,„,,AX) := U^'^tX^^ - j8) / [{'"T)i'n + n+ - j5)]. 

Proof Using (2.2) we write the quantity (2.10) in the form Y,'^^T:7tt,Tcl. Next, consider the 
sequence {wa;;^ = 7r;t;T/'";f;m,n(^)}JL„,+„+i, where r^.„,,„(^) ai^e the numbers given by (2.9), and 
observe that this sequence is increasing in k, with w,n+n+\;T = 1- D 

The upper bounds (when there are at least two) of the residual R,n^n{g), given by (2.10), are 
not comparable. For example, consider the functions = (p^ and g2 = (Pm+n+2 (both belong 
to J^°°{X)), where (pt are the polynomials given by (A.l), and observe that 

u,Kq^X){g['\x))=u,T\n%t2i'^-jS)>0 = u,+,^q'+'{X){g['+'\x)) 

and 

u^'Eq^X){g^^\x)) _ {m + n-T+l){l-{m + n + T)d) ^ 
Ur+iEq'+'{X){gi'+'\x)) ~ im + n-r + 2){l-{m + n + T+l)5) 
since 5 < 0. 

Let X ~ IP(/x;5,j3,7) with 5 < 0. Then X is a linear function of a normal, a gamma or 
a beta random variable, see [2]. The bounds S,nn{g) for the three main cases are included in 
Table 2.1. 



3 Investigating the bounds Sm,n for fixed n 

Next, for n fixed, we investigate the bounds S',,,,,, (g) as m increases. We compare the variance 
bounds Sn,n{g) and 5„, given by (2.8) [for m = n] and (1.2), respectively. Also, we compare 
the new upper variance bounds Sn,i{g) and 5„_ij(g) with the existing Chernoff-type upper 
variance bounds S„ (st^) and S„ (•^^,eaJ^) , respectively; see (1.4) and (1.5). 
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Theorem 3.1 Let X ~ IP(jU;5,j3,7) with 5 < 0. Fix the positive integer n and consider a 
function g G J^'^'^{X), where M can be finite (M > n) or infinite. Then for each m\,m2 such 
that <mi <m2<M the following inequality holds 

\N?ccg{X)-S,nU8)\ > Cm,m2,„(5)|Varg(X)-S^,,„(g)|, (3.1) 

where ^,„. := ^ — — > 1. The equality holds if and only if the function 

g : J{X) -^His a polynomial of degree at most n + m\. 

Proof Consider the positive sequence {C/t;mi,m2,n(5) = rk;mun{^) / n-rmA^)] k>mo+n'' where 
rk;m,n{X) are given by (2.9). This sequence is decreasing in k. Specifically, 

Moreover, we observe that r^-fji^ fi {X) > and 

fk;m2,n (^) — all ^ — ?1 + /Mi + 1 , . . . , « + m2. 

Therefore (3.1) follows. 

We write | Varg(X) -5„„,„(g)| - Cmi,m2,«(5)| Var^(X) -S„,2,„(g)| = 'Lk>n+mi ^kcj and we 
observe that dk > for all k. Thus, the equality in (3.1) holds if and only if g is a polynomial 
of degree at most n + mi. □ 

Notice that if 5 < then U%t"'\^ - jS) /itJitT^C^ " ^'^) > ^ ^^^^ " ^^'^ < "^2- 
Therefore, Cmi,m2,«(5) > Cmi,m2,«(0) = (m2 + n)„/(mi +«)„, since 5 < 0. 

Remark 3.1 Assume the conditions of Theorem 3.1. (a) In view of Remark 2.1(a), the bounds 
{Sm,n{g)}m=o ^^c Same type. From (3.1) it is follows that the bound S,n2.n{g) is better 

than the bound S,„^,n{g). Now, writing n = 2r (when n is even) or n = 2r+ 1 (when n is odd) 
we observe that 

SQ.2r{g)<Sl,2r{g) <■< Varg(X) < ••• <5i,2.+ lte) <5o,2r+lte). 

(b) For the case M = oo, from (2.1), (2.8) and (a) it follows that 

(g) / Varg(X) or Sin.n 
[when n is even] [when n is odd] 

Now, we compare the existing variance bounds 5„( = So.n{g)) with the best proposed 
bound shown in this article requiring the same conditions on g, i.e., with the bound S„„{g); 
see Remark 2.1(b). 

Corollary 3.1 The variance bounds S„ji{g) and Sn, given by (2.8) (for m = n) and (1.2) re- 
spectively, are of the same type and require the same conditions on g. Moreover, the new bound 
Sn.nig) i^ better than the existing bound Sn- Specifically, 

\Y^g{X)-Sn\ > f^;)^|^|Var-,(X)-5„,(,)|. 



The equality holds only in the trivial cases when Varg(X) = S„ji{g) = Sn, i-e., the function 

g : /(X) -^His a polynomial of degree at most n. 

Note that, since 5<0, C:)U%-\^ - j5) /U%'o{l - jd) > f;). 
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The quantities 5„ (^tr) ^^id Sn,i (g) are upper bounds for Varg(X). Both bounds are equal to 
Varg(X) if and only if the function g is a polynomial of degree at most n+ I. The quantities 
S„ (^^eak) ^rid Sn-i,i{g) are upper bounds for Varg(X). Both bounds are equal to Var^(X) if 
and only if the function g is a polynomial of degree at most n. Thus, it is reasonable to compare 
these bounds. 

Theorem 3.2 For n = 1,2, .. . and any suitable junction g we have that: 

(a) S„^i{g) < S„ (str), where the equality holds when n = \ or n> \ and g is a polynomial of 
degree at most n+\. 

(b) 5„_i J {g) < 5,, (weak)' where the equality holds when n = \ or n> \ and g is a polynomial 
of degree at most n. 

Proof (a) From (2.8) and (1.4), via (2.2) and (2.3), we have that 

'JnXstr) •^n,\\,g) - Lk>n+\ {^ir+\){^Y^n8)\ nl\%Xi(^-j^) / 

where each coefficient of c^, ^ > « + 1, is zero when n = \ and is positive when n> \. 
(b) Similarly, from (2.8) and (1.5), via (2.2) and (2.3), it follows that 

'^«,(weak) -^i-LlUj - Lk>n „[l-(„-l)5] 1^ n-";-(l-y5) / ^ ' 

where each coefficient of c^, k> n, is zero when n = I and is positive when n> I. □ 
Remark 3.2 For each ?i = 2, 3, . . . it follows that: 

(a) The bound S„,i (g) is better than the bound 5,, (str); notice that the bound Sn,i (g) requires a 
milder finiteness condition, g G Jif"'^ (X), compared to S„ (st^) which requires that g G Jif"'"(X) 

cjr"'i(x). 

(b) The bound Sn-i,i (g) is better than the bound S'^ (^eak) ; in (a), the bound Sn-i,-i (g) requires 
a weaker finiteness condition, i.e. g G J^"^^-^{X), rather than (^eak)' i-C- g ^ '^"'"{X) C 

jr"-i'i(x). 

Final Conclusion The variance bounds given by Theorem 2.1, for appropriate choices of n 
and m, either provide existing univariate variance bounds or improvements. Our bounds cover 
all usual cases, namely: 

- Chernoff-type, see [13, 3, 8, 5, 14, 17, 2]; 

- Poincare-type, see [15, 11, 10, 1]; 

- Bessel-type, see [4, 1]. 

Notice that no further conditions on the function g are imposed; instead the new bounds require 
the same or weaker conditions, see Remarks 2.1, 2.2 and 3.2, Theorem 3.2, and Corollary 
3. 1. Therefore, the new proposed variance bounds outweigh all the existing variance bounds 
presented in the bibliography. 
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Appendix 



A Some useful properties of the Integrated Pearson Family 

Consider a random variable X with density / --^ IP(/i;;jf) = IP(/x; 5, j8, 7). 

Let a e ( 1 , +°°) . Then £ |X | ° < 00 if and only if5<l/(fl— 1). Notice that X has finite moments of any order 
if and only if 5 < 0; see [2, Corollary A.2]. 

The support of X is the interval J{X) = (a, ft)) and the density / e C°°{a, co), see [2]. 

If E|Xp'+' < 00, / e M* = IN X {0} (that is, 5 < 1/2/), then the random variable X; with density function 
fi{x) = q'{x)f{x)/Eq'(X) follows IP(/i,;^,) distribution with /i,- = {/i + - 2i5) and qi = q/{l - 2/5); see 

[2, Theorem D.2]. Note that if E|Xp' < 00 and E|X|^'+' = 00 then the function fj is a probability density function; 
however, V,\Xi\ = 00 so Xj does not belong to the Integrated Pearson Family. 

If E|A'|^'' < 00, A' e M* (that is, 5 < 1 /{2N — I)), then the quadratic q generates the orthogonal polynomials 
through the Rodrigues-type formula, [16], 

Pk{x) = ^j^-^[q'{x)f{x)], xeJ(X), k = 0,l,...,N. 

Afendras et al., [1], showed an extended Stein-type identity of order n. This identity takes the form lEjPi^{X)g{X) = 
E/(J>(:)^W(X), provided that E/(X)|^(*)(X)| < <x>. Also, EPk{X)P,„{X) = 8t,,„k\Eq^ {X)UfS^li{i - j5); thus, 
the system of polynomials {<Pj:}^q is orthonormal, with respect to density /, where 

(pk{x)=Pk{x)(^k\E^{X) n (l-jS)) ■ (A.1) 

j=k-\ 

Hence, 

/ 2<r-2 \-l/2 

Eq''{X)s^''Hx)=(k\TEq''{X) II {I- jS)) E(pk{X)g{X). (A.2) 

j=k-l 

Moreover, the system of polynomials {(pk+i}^-Q [where (pk are the polynomials given by (A.I) and cpj^'^- is the 
j-th derivative of cpk+i] is orthogonal with respect to density /,-. Specifically, if (pkj are the orthonormal polyno- 
mials corresponding to the density /,■ then <pf'^,(x) = vj^'^ q>kj{x), where v^'' = vl.'\x) := [{k + i)i[Yl'jt:k'+i^ii^ ~ 
j8)) /Eg'(X)] see [2, Corollary D.4]. 

Lemma A.l (see [2]) Let a random variable X ~ IP(/x,g) and consider the strictly positive integers n and N 
such that n<N and E\X\^'^ < <x>. 

Ifge M'"-" {X) then g<'' e (Xi) for each / = 0, 1 1 . (A.3) 

E<P*./ (X/)g(''' (Xi) = v*" (X)-E(pk+i {X)g{X) for each { [= {^^^f;;;;;^ _ i (A.4) 

If the parameter 5 of 5 is non-positive then the moment generating function of X is finite in a neighborhood 
of zero; thus the system of polynomials {<Pt}^^Q forms an orthonomal basis of L-(1R,X) and the Parseval identity 
holds, see [2]. Notice that for each / e 1^ the parameter 5/ = is non-positive too. Thus, the system of 

polynomials {(Pkj}'k^^Q is an orthonomal basis of L^(]R,X,) and the Parseval identity holds too. 

B The solution of the system (2.5) 

Consider the determinants dm^n — ^^^i^uun) ^md = det(A 

i;m,n)t i — lj2, ...,n, where the matrix A,;,„„ is 

formed from A,„_„ by replacing column / with the vector 1„. 

For each t = 1,2, ...,n define the matrix S„,.„(f) £ ]r("-'+i)x("-'+i) [„, 

,n are fixed] which has elements 
br,c;m,n(t) = (m + r— l)c- 1 n','^m+r+,Ii ( 1 — jS), where empty products are treated as one. Observe that 
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rf,„,„ = (m + «)„(n"5r'(l-i^))det(fim,«(l)) and 

det (S,„,„ (f ))=(„- f) ! ( ( 1 - [2 j + (f - 1 )] 5)) det (S,„,„ (f + 1 )) , f = 1 , 2, . . . , n - 1. 
Thus, it follows that 

d„^,n = {m+n)„ [xiyj^p] [n}L+r'(i -i^)] n;'r/(n';L+',7'i(i - [2./+ - 1)]^)) ^ 0. 

Now, for each t = 1,2, ... ,« define the matrix Bi^,„j^{t) e ]l^(n^'+l)><("^'+i) [i^m,n are fixed integers] which 
has (r,c)-element br,c-i,m,n(t) = ('" + '')c-i n7=m+r+f-2(l ^ -''^)' "^^^^ ^= \,2,...,i~t, and br,c;i,m,n(t) = (™ + 
'■)c I\l]l!iXrX't-2 ( 1 - ./^) ' when c = - r + 1 , - f + 2, . . . , n - / + 1 . Observe that 
rf/;«,,« = (-l)'"'det(B,v„,„(l)), 

det (5,-,„,„ (f )) = ^pij^ ( n7=+";'i ( 1 - [2; + (f - 1)] 5)) det (f + 1)) , f = 1 , 2, . . . , / - 1 , and 

det (5,-,„,„(0) = (n - . + 1 ) ! (n^j^;; ( 1 - [2.7 + (/ - 1)] 5)) [m + n - z)„_,- [Y{%:,+]{ 1 - det (S„,,„(z + 1 )) . 
Thus, it follows that 

rf,;«,,„ = (-ir'^^^^±HiF[n;u;!] [n7i:;;/(i-75)]nr=/(n:;L+';;'i(i-[2i+(^-i)]5))- 

Therefore, according to Cramer's rule, (2.6) follows. 



C A useful hypergeometric series 

Lemma C.l Let m,n,k (E IN and 5 <0. Then the following hypergeometric series holds: 



ij {m+n)i n;'L+;r'(i (™+")«n"L+r'(i 



(C.l) 



Proof For the case 5 = 0, write (C.l) as 1)'(") (in+'n) ~ ^"{m'+n)^" observe that this follows from 

Vandrermonde's formula L|'=o(^l)' (") {x+y)- ~ (t+y) ' replacing x with k and y with m + n — k, see [6, p. 125]. 
For the case 5 < 0, write (C.l) as £«^„(-l)' ^ . This follows from Dougall's 

Identity, 1^0 „l,4il^ff;+.), = "^'ng substitution a^k,p^{l/S + l-k), {-,n- 

n — 1) and s n> n, see [9, eq. (2)]. □ 
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